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! Abstract 

o 

(N 

' The quantum-classical limits for quantum tomograms are studied and compared 

■ with the corresponding classical tomograms, using two different definitions for the 
i limit. One is the Planck limit where in all ^—dependent physical observables, 

and the other is the Ehrenfest limit where ^ — > while keeping constant the mean 
value of the energy. The Ehrenfest limit of eigenstate tomograms for a particle in 
a box and a harmonic oscillator is shown to agree with the corresponding clas- 
I sical tomograms of phase-space distributions, after a time averaging. The Planck 

■ limit of superposition state tomograms of the harmonic oscillator demonstrates the 
^ . decreasing contribution of interference terms as ^ — > 0. 
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1 Introduction 



The quantum-classical transition has been the subject of numerous studies, 
[1,2,3] for the classical limit of bi-Hamiltonian systems see e.g. [4]. In the 
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standard formulation of quantum mechanics the states of a system are as- 
sociated with wave functions[5] or density operators, [6,7] and the WKB pro- 
cedure can be employed to take the classical limit. Further insight in the 
connection to classical mechanics can be gained from other formulations of 
quantum mechanics, [8] e.g., the Wigner-Moyal[9,10] and Feynman path inte- 
gral formulations. [11] In the Moyal approach, the Wigner function satisfies 
the quantum evolution equation which is similar to the kinetic equation for 
classical probability distributions in phase-space. In the limit h ^ the Moyal 
equation yields the Liouville equation of classical statistical mechanics. The 
classical action plays an important role in the Feynman path integral represen- 
tation for the Green functions (propagators) of the Schroedinger equation, and 
the only contributions that appear in the complex probability in the H ^ 
limit are associated with classical paths. Also the classical action satisfying the 
classical Hamilton- Jacobi equation appears in the wave function after using 
the WKB decomposition. 

Despite all the utility of these well known approaches, because objects such 
as wave functions, Wigner functions or other representations of density oper- 
ators(e.g., Sudarshan-Glauber singular P-quasidistribution[12,13] and Husimi- 
Kano Q-quasidistribution[14,15]) differ in an essential way from classical phase- 
space probability distributions, it is not easy to make a precise comparison 
between the two systems. This can be rectified by going to the tomographic 
approach where both the classical and quantum theory offer a discription in 
terms of tomographic probabilities. [16,17,18] The quantum tomogram realizes 
a specific version of star-product quantization [19] and allows for a comparison 
with the classical tomogram, both being written on the same domain. The 
advantage of a tomogram is connected with its property to be a standard pos- 
itive probability distribution function describing the quantum state and this 
property was found useful to study the entanglement criterion for continuous 
variables, [20] as well as to formulate a separability condition for spin degrees 
of freedom. [21] 

The aim of this work is to consider the quantum-classical transition within 
the framework of the tomographic approach. We shall mainly be concerned 
with examples and postpone more general considerations to a future paper. In 
particular, we shall consider the tomographic approach for the description of 
quantum states for a particle in a box, as well as the quantum tomograms of 
harmonic oscillator coherent states and stationary states. We study two kinds 
of quantum-classical transitions, the Planck limit and Ehrenfest limit of the 
quantum tomograms. In Planck limit we obtain the quantum tomograms of 
harmonic oscillator coherent states and stationary states for ?i ^ 0. In the 
Erhenfest limit we obtain the limit of the quantum tomogram of harmonic os- 
cillator stationary states and of the states of a particle in a box for /i — > with 
fixed energy. We will show that the Ehrenfest limit provides the time-average 
of classical motion (both for the harmonic oscillator and the particle in a box) . 
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Also we study the behaviour of tomograms for superposition states of the har- 
monic oscillator in the Planck limit. Namely, we show that the tomogram of 
the superposition of two arbitrary eigenstates (the Schroedinger cat[22]), goes 
to a mixture of classical tomograms of these states. The interference contri- 
bution to the tomogram goes to zero in the Planck limit. Another type of 
Schroedinger cat system (consisting of even and odd coherent states[23]) will 
be considered in both the Planck and Ehrenfest limits. 

The article is organized as follows: in section 2 we review known properties 
of symplectic tomograms for classical and quantum systems. In section 3 we 
introduce the quantum-classical transition. In subsection 3.1 the Planck limits 
are computed for few systems. In subsection 3.2 we study the Ehrenfest clas- 
sical limit for the stationary states of the harmonic oscillator and the infinite 
square well potential. Concluding remarks and perspectives are given in the 
final section. Some relevant formulae for generalized functions and their limits 
are considered in the appendix. 



2 Symplectic tomograms in classical and quantum mechanics 

Following references[16, 17,24] we review the tomographic description of par- 
ticle states in classical and quantum mechanics. For simplicity we restrict to 
one-dimensional particle systems. 

2.1 Classical tomograms 

The standard description of classical states with fluctuations is given by a 
non-negative joint probability distribution function /(p, q) on the phase space 
of the particle with one degree of freedom. The function is normalized, i.e. 



The tomogram is the probability distribution function in a rotated and scaled 
reference frame on the phase-space. The classical tomogram is constructed 
from /(p, q) and is a function of a coordinate X e M, which is related to the 
position q and momentum p of a canonical reference frame on phase-space by 



where and u arc real parameters. It can be expressed in terms of a scaling 
parameter s and a rotation parameter 9 : 





X ^ ljLq + up, 



(2) 



II — s cos 9 , y — s ^ sin 9. 



(3) 
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For fixed fi and u one then gets a line on the commutative plane (g, p) with an 
orientation 6 from the position axis. Thus the physical meaning of the variable 
X is that it is the "position" of the particle measured in the reference frame 
of the phase-space whose axes are rotated by an angle 9 with respect to the 
old reference frame, after preliminary canonical scaling of the initial position 
q ^ sq and momentum p — > s~^p. The coordinate X of equation (2) together 
with 

Y = -s^vq + s^'^jip (4) 

provides a canonical transformation preserving the symplectic form in the 
phase-space. For that reason the classical tomogram is called "symplectic" . 

The tomogram of the classical statistical density /(p, q) is defined by the 
Radon transform[25] 

/X, i/) = j f{p, q)S{X -i^q- iyp)dpdq. (5) 

It can be written in the form 

WiX,i,,u) = {6iX-f^q-up))^ (6) 

where the average is done using the probability distribution f[p, q) in the 
phase space. Its inverse is 

f{p, q) = 1^-^ J >V(X, n, v) exp [i{X - fxq- up)] dXdndv. (7) 



Since the classical probability distribution f{p,q) is normalized, the classical 
tomogram is also normalized for any values of the parameters /j, and u, i.e. 

Jw{X,ij,iy)dX^l. (8) 

For II — 1 , 1/ — the tomogram provides the marginal distribution of the 
position {X — q) 

W{X,1,Q) = I f{p,q)dp. (9) 

For = , z/ = 1 the tomogram provides the marginal probability distribution 
of the momentum [X = p) 

>V(X,0,1)= j f{p,q)dq. (10) 

The expression (5) is written for the case of stationary distribution function in 
phase space. But it is easily generalized to time-dependent functions /(p, g; t) 

W{X, II, i^;t)^l f{p, q; t)5iX -iiq- vp)dpdq. (11) 
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For distribution functions associated with particle motion along a given tra- 
jectory in phase space q{t),p{t) 

f{p,q;t)^S{p-p{t))S{q-q{t)) (12) 

the tomogram given by Eq.(ll) reads 

Wf{X, ly; t) = S{X - f,q{t) - vp{t)). (13) 

For example, classical free motion with initial position go and initial momen- 
tum po, i-e. q{t) — qo +Pot] pif) = Po, is described by the tomogram 

Wf{X, IX, p;t) = S{X- ix{qo + Pot) - i^Po). (14) 



Below we will compare tomograms associated with motion along classical tra- 
jectories with classical limits of their corresponding quantum tomogram. When 
the trajectory is periodic it is possible to introduce the time-averaged classical 
tomogram 

1 r^. 



(W) (X, /X, z/) = - / W(X, /i, t)dt. (15) 
I Jo 

The time-averaged tomogram corresponds to the time-averaged distribution 
function on phase-space 

{f)ip,q)-^£ f{p,q;t)dt. (16) 
2.2 Quantum tomograms 



Given a density matrix p[x, x'), one constructs the corresponding Wigner func- 
tion W as: 

W{p, q)^ I Piq+'^^q- |) exp(-i^)dii, (17) 

and the Radon transform of Wigner function W is the quantum tomogram W 
of p, which is a positive function[18] of three variables, 

W{X, ^) = ^ / W{p, q)5iX -n- up)^ (18) 



1 r u u 
Equation (18) can be written in a form similar to Eq.(6), 



.X - nq 
—t — u 



dqdu (19) 



W{X,ii,v)^{5{X-m-vp)). (20) 

The difference with Eq. (6) is that here the position and momentum are quan- 
tum operators q and p , i.e. we use non-commutative geometry of phase-space 
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plane, taking into account uncertainty relations. The averaging in Eq.(20) is 
done using density operators p, i.e. 

(21) 



A) := TrpA. 

The inverse formula of Eq.(18) is readily written as: 

W{p, q) = ^ [ W{X, 1^, u) exp [i{X - jiq- up)] dXdjidu. (22) 
27r J 

The density matrix p(x, x') can be obtained from the Wigner function as 



p{x,x') 



2TTh 

In terms of tomogram we get: 



x + x' 
W{p, — - — ) exp 



p(x — x') 



dp. 



(23) 



p{x, x') = ^J W{X, p, exp 



h 



^iX-p^) 



dXdp (24) 



For pure states p{x,x') = ip{x)ip*{x') and W can be computed from Eq.(19): 
yV{X,P,iy) = / i'iq+'^Wiq - '^)e-'^''-^^^ndqdu 



1 



271 h \ u\ 



(25) 



Thus, apart from the pre-factor l/2TTh \ the tomogram W of the wave func- 
tion ip is the square modulus of the tomogram amplitude : 



A^{X,p,iy) : = 
When u — and p one can instead use 



W{X,p,v)^ 



j 'lp{p)e '2n^P^+''n,,P(lp 



2nh \p\ 

where the Fourier transform of the wave function has been introduced 

1 



V'(p) 



i>{y)e '^dy. 



(26) 



(27) 



(28) 



If both p — 1/ — 0, we get W(X, 0, 0) = S{X), after using / p(x, x)dx — 1. 

Tomograms of quantum states and classical statistical density are positive 
functions on the same space and therefore can be compared. Moreover, the 
same dimensional relations hold in both cases: 



[W{X,p,iy)] = [X]-' 



(29) 
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where [.] indicates the units. X and W can be made dimensionless after as- 
suming 

M = M-^ ; M = b]-^- (30) 



The operator 

X = nq + up (31) 

together with its conjugate 

Y ^ -s^uq + s''^HP (32) 

preserves the canonical commutation relations: = [q,p\- The observ- 

able X is a new position operator, i.e. the position after a symplectic (linear 
canonical) transformation in the quantum phase-space {q,p) of the particle. 
The real variable X gives the possible results of a measure of X and runs 
over the spectrum oi X. In this way a description of quantum tomograms is 
recovered in complete analogy with the classical case. So the tomogram is also 
"symplectic" in the quantum case. 

In classical mechanics, the transition from the distribution function of two 
canonically conjugate variables (position q and momentum p) to the distri- 
bution function of one variable (position X) does not play a crucial role due 
to absence of quantum mechanical constraints like the uncertainty relations 
of Heisenberg[26] and Robertson-Schroedinger[l,27,28,29,30]. On the contrary, 
the use of tomograms in quantum mechanics provides the possibility of describ- 
ing a quantum state with a probability distribution which does not violate the 
uncertainty relations. Comparing the Eqs.(6) and (20) we see that the classical 
limit will be found starting from a non-commutative plane {q,p) and going to 
a commutative one {q,p)- This means that in the limit process we take into 
full account the behaviour of the quantum uncertainty relations. 



3 The quantum-classical transition 

The classical limit of quantum states can be obtained using different proce- 
dures. One procedure is to take the limit ^ of certain quantum mechan- 
ical expressions. The wave function defined on the configuration space is not 
suitable for this purpose, since there is no physical interpretation of the re- 
sulting limit. It is more useful to associate an "observable" with any state 
vector, namely the associated rank-one projector. Following the Wcyl-Wigner 
approach observables can be mapped to functions on phase-space and finally 
to tomograms via Radon transform. The classical limit h ^ then can be 
done on tomograms in two different ways. 
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(1) At a kinematic level (Planck limit). Here no reference is made to any 

specific dynamics. 

(2) At a dynamical level (Ehrenfest limit). Here the limit ^ — > on tomo- 
grams is performed keeping constant the mean value of a given observable 
(usually the Hamiltonian) . 



3.1 Planck limits 



Wc can anticipate some general features of the Planck limit for the tomogram 
of a quantum state. The quantum fluctuations of position and momentum, 
depending on Planck's constant, are washed out when — > and in general a 
Dirac 5{X) results, corresponding to the tomogram of the classical probability 
distribution 5{p)5{q). The fact that that Dirac delta function is centered at the 
origin has no particular physical significance, as it is only due to the arbitrary 
choice of the origin of the affine phase space. 

The Planck limit can be predicted from a scaling argument for a large class 
of wave functions. Assume the dependence on h in the wave function ■0 is (7 
a real number) 

^(x) = him{K'x). (33) 

This assures that 

j \i^{x)\^dx^l (34) 
independent of the value of h. Computing the tomogram in this case we get 



27r/i \v\ 
1 



27r|i/| ^T+i 



'^*dt 



(35) 



When 7 = —1/2 it gives 

1 



27r \v\ Vh 



(36) 



Since in general J yV{X, fi,i')dX = 1, the Dirac delta theorem of the Ap- 
pendix, with n — 1/ Vh, applies and yields the Planck limit of the tomogram 
as 

limyV{X,pi,u)^5{X). (37) 

ft— »u 

This is the case of the harmonic oscillator eigenstates and their superpositions. 



When 7 = for the wave function in position space, its Fourier transform ijj of 
Eq.(28) scales with 7 = —1/2 and using Eq.(27) we recover the same Planck 
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limit of Eq.(37). In general this happens when a quantum particle of mass m 
is confined in a potential with a characteristic length L. then Vq = /2mL'^ is 
the scale factor of the potential energy so that the eigenvalues are proportional 
to Vq while the eigenstates do not depend on h. This is the case of the Poeschl- 
Teller potential[31] and of its hmit, the infinite square well. 

We finally observe that a more general scaling law could be 

'il;{x)^K^^{W{x-Xo)). (38) 

In that case, when 7 = —1/2, the Planck limit of the tomogram is 5{X — iixq) 
and corresponds to the classical distribution 5{p)5{q — xq). Nevertheless, by 
using the shifted position operator q — xq, the Planck limit turns out to be 
again S{X). 

In the following we study the Planck limit for examples involving Hermite 
eigenfunctions. We start by reviewing the Hermite tomograms, whose expres- 
sions will also be useful in evaluating the Ehrenfest limit. 



3.1.1 Harmonic oscillator 



Recall that the eigenstates of a quantum harmonic oscillator of mass m and 
frequency uj, normalized with respect to the Lebesgue measure dy, are 



^n{sl^y) = (x/^2"n!)-^^-i/„(^-|/)exp 



(39) 



with w = muj and n = 0, 1, 2, ... Using Eq.(25) the tomogram yVo{X, ji, v) for 
the ground state can be computed as 



1 



exp 



-X' 



Similarly, the tomogram >V„ for the n-th excited state is 



X 



/i(roV2 + /x2) 



X 



(40) 



(41) 



A suitable way to obtain these results is to evaluate a generating function J 
for the tomogram amplitudes of Hermite functions 



s s r w 

n=o V n! „=o yn\J V ^ 



exp 



— 1/ — t—y 
. 2ni' hi' . 



dy 



w 1 2T^hv 




—\l exp 



C 2 . 

2C ' 



l2wXs 



__ _ X 

' ^ 2huC* 



(42) 
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where ( = zuiy+ijj, and (* is its complex conjugate. It is possible to reconstruct 
a generating function of Hermite functions Defining 



(43) 



Equation (42) becomes 



c 



J(t) — v2T^hu\l — exp 



oo 



n=o vn! 



(44) 



Eventually, the Hermite tomogram amplitudes can be written 



1 



(0) = \/2T:hu exp 



i-A-'4\ MQ), (45) 



which yield Hermite tomograms W„ of Eq.(41) for (/?„. Observing that is 
normalized to unity, the Planck limit of Wn, evaluated by means of the Dirac 
delta theorem (see Appendix), is 



]imWn{X,fx,iy) = 6{X) 
h—*o 



(46) 



as expected. 



3.1.2 Superposition 

The case of a superposition of two harmonic oscillator eigenstates 



(47) 



shows in the simplest manner the vanishing of quantum interference in the 
classical Planck limit. Writing the tomogram in terms of Hermite amplitudes 
of Eq.(45) 



1 



2TTh \u\ 



1 

72 



i^n + ^m) 



A A* 
Reli!^. (48) 



Dropping a unimodular factor independent of h, the quantum interference 
term reads 

^^,^Vn{Q)^m{Q). (49) 

inh \v\ 

Due to the orthogonality of harmonic oscillator eigenstates, this term vanishes 
when — > as is shown in the Appendix. Moreover, it is possible to show that 
the quantum interference term vanishes like \/h. Again, by using Eq.(46), the 
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classical Planck limit of superposition state tomogram is 



\imW^{X,ii,iy) = S{X). 
h—^o 



(50) 



3.1.3 Coherent states 

We next consider the tomogram of coherent states. Recalling that the coherent 
state I a) labeled by the complex number a is defined by 



n=o vn! 



(51) 



the tomogram amplitude of the coherent state can be expressed through 
Hermite amplitudes of Eq.(45) 



— e 



-|"l'/2 \- 



n=o V n! 



-H- 



/V(a) 



(52) 



From the expression for the generating function of Hermite amplitudes J in 
Eq.(42), the previous equation becomes 



A„ = — — : — e 'aw^*^ e 2 exp 



TT^V C 

SO that, eventually, the tomogram Wq(X, i/) — AaA'^/27rh\u\ is given by 

(^^/wX — ij,\/2hRea — wv^/2h\ma 



(53) 



exp 



(54) 



After shifting the integration variable 

j Wa{X, u)dX = J^exp [-e] dC = 1, 
and by the Dirac delta theorem the Planck limit of the tomogram is 



}Va{X,l,,u)-^d{X). 

n— ♦0 



(55) 



(56) 



3.1.4 Schroedinger cat states 

Even and odd coherent states (Schroedinger cat states) are defined[23] as: 

\a±)^N±{\a)±\-a)) (57) 
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where 



I a) = exp 



N± = , ^ (58) 
/2(l±e-2|"P) 



2 2 

The tomograms yVa±{X, /x, i/) = A^^ \Aa ± ^-a|^ /27r/i contain three terms 

Nl {Wa{X, u) + >V_«(X, /X, u) ± 7(X, I/)} . (59) 
The interference term reads 

I(X,^i,u)^2Re^^ (60) 

ZTTft |Z/| 

with tomogram amphtude Aa given by Eq.(53). As shown in the Remark of 
Appendix, it results 



/ 



^^^dX = {-a\a) = e-2H' (61) 



and by Dirac delta theorem the Planck limit of the interference term /(X, /i, u) 
is 2e~2'°' S{X). Then the Planck limit of the Schroedinger cat tomogram is 

\imWaJX,ii,u) = Nl(25(X) ± 2e-^l"l'5(X)) = 5(X). (62) 

/i— >0 

We observe that the Planck limit S{X) is gained thanks to the contribution 
of a non-zero interference term between non-orthogonal states. 



3.2 Ehrenfest classical limits 



The Ehrenfest classical limit takes place at the dynamical level, i.e. it requires 
a specific dynamical system and the limit is performed while keeping constant 
the value of some selected observable, usually the energy. [1] In this section 
some Ehrenfest classical limits are evaluated with the constraint of constant 
energy. 

We observe that after calculating the Ehrenfest limit, it is possible to consider 
the hmit when the energy goes to zero. This is not equivalent to the Planck 
limit because in the Ehrenfest case, if there is a potential with a minimum, the 

particle goes to its equilibrium position go when the energy vanishes. Then the 
expected result is 6{X — fiqo), rather than S{X), agreeing with the classical 
tomogram of a rest state 6{p)6{X — qq), as the next example of the coherent 
state shows. 
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3.2.1 Schroedinger cats 



The Ehrenfest classical limit of the coherent state tomogram Wa is readily 
evaluated from Eq.(54). Put zu — 1 and impose the constraint on the energy 
by taking constant the mean value of the harmonic oscillator Hamiltonian 
h{a'^a + 1/2) on the coherent state \a). The constraint is satisfied by choosing 

aVh = const =^ ha* a = h (a^aj = const . (63) 



Now, bearing in mind that the mean value of the annihilation operator a on 
the coherent state \a) is 

(a)^ = - — , = , = a , 64 



the Ehrenfest classical limit — > , a — > oo with ay2h — const — Qa + iPa 
turns out to be, just by inspection of Eq.(54) (with w — 1), the displaced 
delta function 

S{X - fxqa- UPa) (65) 
which coincides with the tomogram of the classical distribution function 

f{p,Q) = S{p-Pa)S{q-qa)- (66) 



As far as the Schroedinger cat states of Eq. (59) are concerned, the interference 
term I{X, ji, v) vanishes in the Ehrenfest hmit, due to the presence of a rapidly 
oscillating factor (/i — > 0) 



2v2^-X(z/Rea — ulma) 2X (lyqa — upa) /^„s 

while the normalizations N"^ go to the same constant value 

^± = 2(l±e-2|"P) ™ 2 ^^^^ 
and the tomograms become the same sum of delta functions 

Woo{X, H, u) = ^[S{X - liqa - VPa) + 5(X + /Xg^ + VPa)]- (69) 



3.2.2 Infinite square well 

We consider a particle of mass m = 1 in an infinite square well < g < L. To 
deal with wave functions of compact support it is convenient to regard them as 
elements of L2(M) suitably projected by means of the projector X[o,l] , which 
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is the standard characteristic function of the set [0, L] . The eigenf unctions and 
the eigenvalues are 



^n{q) = ^sin(^g)x[o,L](g) ; E., 



n-^i^f, n=l,2,... . (70) 



Fixing the energy to be one for convenience, the Ehrenfest classical limit 
amounts to 

- - (71) 



n — > oo. 



In general, for u 0, 



rnr 



/ sm(— |/)e 2fi-^ ^fi-^dy 
Jo L 



where the amplitudes are respectively given by 



n 


A_-A+ 


2 


AL^V2 







(72) 



J exp[inF:^{y)]dy J explain— 



X 



2a/2, 



-y - 



V 



V2, 



V 



dy. (73) 



These amplitudes can be computed exactly in terms of Error function. How- 
ever, in view of the large n limit, it is more useful to evaluate these integrals 
in the stationary phase approximation. It results 



A^ ~ 



2L\v\ v/2 
11 



■ exp < I 



nF^(g?)+sign(^)^]} Xio^m- (74) 



where the stationary points are 



(75) 



Eventually, the asymptotic behaviour of W(X, /x, v) is 

X[o,L](Q7) + X[o,L](Q+) - 2x[o,L](g;)X[o,L](Q+) cosn[F_(Q7) - F+(Q+)] 

2|/x|L 

and evaluating the Ehrenfest classical limit n — > oo in the sense of a distribu- 
tion: 

w^{x,ii,v) = 2^ {x[o,L](g;) + x[o,L](g+)} . (76) 

This result has to be compared with its classical analogue. In that case, bearing 
in mind that m — 1 and p — ±-\/2!B = ±-\/2, we have 



where q{t) is a zigzag line of height L and half-period Lj \f2^ 

q{t) = V2tX[o,L/V2]ii) -^{t- V2L)x[L/^2,V-2L]it), 



(77) 



(78) 
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while 

Pit) = V2X^0,L/V2]it) - ^X[L/V2,V2L]{t)- (79) 

To make a comparison with the Ehrenfest hmit we have to average the classical 
time dependent distribution tomogram over a period: 

1 /-V^L r 

W{X, II, u) = dt J fip, q; t)5{X -iiq- yp)dpdq 

I ( i-y/2L/2 rV2L \ 

= 71^;^ ( io " ^^^^^'^^ + W ^^^^ " ^^^^^'^V ' 

where are again given by Eq.(75). Eventually, 

>v(x,//,i/) = {x[o,L](g:) +x[o,L](g+)} = >Voo(^,//,i^). (si) 

This general result yields the marginal probability distribution of g, in the 
limit z/ — > 0, yU. = 1 

{x[o,L](g;) + x[o,L](g.+)} = \xyo,Lm. m 

As expected, this result shows that the particle position is likely distributed 
in the box. 

The limit — > 1,// = which yields the marginal probability distribution of 
p is more involved. However, observing that 

/ ^X[o,L](g?)ciX = l (83) 

and remembering the explicit expression of QJ, by Dirac delta theorem it 
results 

JS-^o^^ {x[o,L](g;) + X[o,L]iQt)} = l^iX -V2) + h{X + V2) (84) 

As expected, the probability distribution of momentum is likely concentrated 
on the two allowed classical values of p. 



3.2.3 Harmonic oscillator 

We assume m = l,a; = 1 and evaluate the Ehrenfest classical limit of the 
Hermite tomogram Wn of Eq.(41) when 

h^O,En^ const = 1 =^ (n+^ 1. (85) 
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In view of the large n limit, we may assume n + | ~ n, so that h — 1/n. 
Besides, due to the spherical symmetry on the phase space, it is enough to 
consider the case = 1, = 0. Then the tomogram of the harmonic oscillator 
eigenstate ipn is written as 



>V„(X,1,0) 



n e 



lX2 



(86) 



Note that the tomogram is an even function oi X, so we only need study 
positive values of X. We start by expressing Hermite polynomials in terms 
of parabolic cylinder functions U (see [32]) 



(87) 



The main asymptotic formula is 

, 1 a 



T 



U{a,x) 2-VW2r(_ _ ^)^^^Ai{T), (x > 0,a ^ -oo) , 



where the Airy function Ai and the Euler gamma function F appear, with: 

2/3 



X 



2^ \a 



; r=(4|a|) 



2/3 



(89) 



and correspondingly: 



i[arccose-ev/r^], < 1) 
j[eVF^-cosh-i^], (e>l) 



(90) 



In our case, the asymptotic behaviour of Wn{X, 1,0) is 



e-1 



At\r) (91) 



Now, the numerical pre-factor can be readily estimated as: 



7rn! ^ 2 ^ 



n „^r^M ~_ ^^J- = x/2 (92) 



r(f + 1) n-oo 



n 



Using definitions in equations (89), (90) we have 



--X 



X 





3 







3/2 



16 



and the asymptotic behaviour of the Airy function [32]: 

Ai(-|r|) ~7r-i/2|^|-i/4gij^ 
Ai{\r\) ~ (47r)-V2|r|-^/^exp 
we obtain the Ehrenfest classical distribution hmit of W„(X, 1,0): 



2 I |3/2 ^ 

3 I ' I "T 4 

_2 U|3/2 
3 I' I 



(93) 



V2 



sm 



(n + |)(arccose-ev^r^) + - 



k (e<l) 

(e>i) 



4 exp [-2(n + I) (^VF^ - cosh-^ ^ 
with ^ = X/ \/2. Eventually, restoring the negative values of X, we have 

V2 



>Voo(X,l,0) 



2tt, 



X[-Ml(^)- 



(94) 



The classical tomogram yV{X, 1,0) is the time average of the Radon transform 
oi f{p,q;t) : 

i £dt I 6{p-p{t))6{q - qitMX - q)dpdq = ^ 5{X - q{t))dt (95) 

where T — 2t: and q{t) = go cost + posint. Due to the rotational invariance 
of the harmonic oscillator Hamiltonian H — ip"^ + q^)/^ we may assume 

q(t) = qo cos t; p{t) = -qo sin ^ 2 - gg ^ =^ = ±V^ (96) 

so that the classical tomogram W(X, 1,0) is 



27r 



go sin(arccos ^) 



X 1 



X 



1 - tt)^ 



X[-i.il(^) (97) 



which coincides with the quantum Ehrenfest limit of Eq.(94). Furthermore, 



dq{t) = —go sin t dt 



dt 



dX 



|gosint| 



(98) 



and dt/dX is the inverse of the classical velocity V, so finally Eq.(97) can be 
written as: 

2 1 ,X 



W(X,1,0) = 



T 



X[-i,i](^)- 



(99) 



As expected, the probability distribution of the position depends on the inverse 
of the classical velocity V in that point and diverges at the turning points 
X/V2^±l. 
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4 Conclusions 



To conclude we found that the classical limit of quantum state tomograms 
is suitable for a comparison with classical state tomograms. Prom the Planck 
limit h ^ 0, the tomograms of stationary states of the harmonic oscillator, 
as well as of its coherent states, yield the localized state S{X), agreeing with 
the classical tomogram of a state localized in the phase space. The Ehrenfest 
limit of quantum tomograms, where the energy is fixed while — > 0, gives 
the expected expressions of classical tomograms of classical states. When the 
energy vanishes, the result is a rest state 5{X — hqq) with the particle sitting 
at the minimum of the potential. 

The same results were shown to apply in the case of both the Planck and the 
Ehrenfest limits of the quantum tomogram for a particle in a box. 

We have also found that, for the Schroedinger cat states, the interference 
term of a superposition of two orthogonal states vanishes in the Planck limit 
as Vh. In the Ehrenfest limit, even and odd coherent states yield the expected 
mixture of two classical delta distribution. 

For composite systems, we believe that tomographic description of quantum 
and classical states can be a suitable tool to study the classical limit of en- 
tangled states both in the Planck and in Ehrenfest limit. This aspect will be 
considered in a forthcoming paper. 

Finally we observe that the diverging inverse velocity factor in the Ehrenfest 
limit of the harmonic oscillator results from the presence of turning points of 
the harmonic oscillator motion along the X— axis. These divergences should 
disappear when describing the motion with a complex variable instead of a 
real one. To do that, the holomorphic (i.e. Bargmann-Fock) representation 
is suitable. However, this requires a preliminary discussion for dealing with 
tomograms associated with non-Hermitian operators, which will appear else- 
where. 



5 Appendix 

Many results of the present paper follow from a Theorem, valid under quite 
general assumptions, that here we recall: 

Dirac delta Theorem Let / be a summable function on the real line such 
that 




(100) 
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Then nf{n{x — x')) — > NS{x — x') when n — > oo. 

Proof Let ip be any test function. Then by Lebesgue Theorem 

hm n j f{n{x — x'))ip{x)dx — hm J f{s)(p{ — \- x')ds 

= J \im^f{s)ip{^ + x')ds = (p{x') J f{s)ds = N(p{x'). (101) 

RemEirk We can evaluate the Planck limit of the product of tomogram am- 
plitudes A^A^,when the state vectors {y\ip) and have scaling exponent 
7 = —1/2. In general, for any 7, it is 



27rh lul 



= J ilj{y)(j)*{u)e'^^y^-''''^5(y - u)dydu = (102) 
In the case 7 = — 1/2, scaling the integration variables yields 

= / ^il{y)^*{u)e^^^(y'--\-'±^y--^^y^. (103) 
27rn.|i/| J 2Trvh\i'\ 

So, as a consequence of the previous Theorem (withn = 1/y/h), it results 

For the tomogram amphtudes An{X, /i, u) of the eigenstates ipn of a quantum 
harmonic oscillator, it results 

hm^^^ = 6n,mS{X). (105) 
where 5n,m is the Kronecker delta. 

When the scahng exponent of {y\ip) , {y\(/)) is 7 = 0, the same result can be 
obtained trough their Fourier transforms which scale with 7 = —1/2. 
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